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MAGNETO-TELLURIC  STUDIES  IN  CENTRAL  ALBERTA 


ABSTRACT 

Magne to- telluric  records  were  taken  at  five  locations 
in  central  Alberta.  The  data  were  analyzed  by  the  method 
outlined  by  Blackman  and  Tukey  (1958)*  The  data  processing 
was  done  entirely  by  electronic  means.  An  interpretation  was 
carried  out  by  employing  Cagniard’s  theory  in  conjunction 
with  information  about  the  substructure  obtained  from  other 
geophysical  and  geological  methods. 

For  the  sedimentary  basin,  resistivities  ranging  from 
6  to  50  ohm  meters  were  obtained.  Moreover,  the  results 
indicate  this  basin  to  be  horizontally  isotropic.  For  the 
Precambrian,  resistivities  ranging  from  600  to  10^  ohm 
meters  were  obtained.  The  results  indicate  the  Precambrian 
to  be  horizontally  anisotropic.  No  electrical  discontinuity 
was  observed  at  depths  ranging  from  l\.0  to  5>0  km.  This 
suggests  that  there  is  no  electrical  discontinuity  at  the 
Mohorovicic  discontinuity.  From  seismic  evidence,  the  depth 
to  this  discontinuity  was  determined  to  be  43  Em  at  Irricana. 
However,  the  results  are  suggestive  of  an  electrical 
discontinuity  at  a  depth  of  approximately  80  km. 


ACKNOWLEDGMENTS 


I  wish  to  thank  the  following: 

Dr.  K.  Vozoff  for  his  supervision  and  assistance  with 
various  phases  of  this  project. 

Dr.  G.  D.  Garland  for  his  helpful  suggestions. 

R.  M.  Kills  for  processing  the  field  data. 

M.  D.  Burke  who  designed  and  assembled  the  field 
apparatus . 

P.  Viney  for  assisting  with  the  drawings. 

Me mb e is  of  both  the  Computing  Centre  and  the  Geophysics 
Group  (University  of  Alberta)  who  have  assisted  with  various 
phases  of  this  work. 

The  Alberta  Society  of  Petroleum  Geologists  for 
permission  to  reproduce  a  regional  geological  cross  section 
of  the  western  Canada  sedimentary  basin. 

The  Petroleum  Research  Foundation  of  the  American 
Chemical  Society  and  the  National  Research  Council  for 
supporting  this  project. 


TABLE  OF  CONTENTS 


Page 


Chapter  I:  Introduction .  1 

Chapter  II:  The  Magneto-Telluric  Method  (Cagniard  1953) 

Outline  of  Method . if 

Frequency  Range  of  Interest  .  5 

Uniform  Fields .  6 

Classification  and  Origin  of 

Micropulsations  .  •  8 

Layering  of  Earth . .  .  11 


Displacement  Current  in  a  Conductor.  .  .  lif 

Direction  of  Propagation  of  Magneto- 

Telluric  Wave  in  the  Earth  .  Ill- 


Intrinsic  Impedance  of  the  Eiarth 
Which  is  Treated  as  a  Semi- 
Infinite,  Uniform  Conductor . 16 

Apparent  Resistivity  .  ...  17 

Two-Layer  Earth  Model . 18 

Five-Layer  Model  .  21 

Interpretation  of  Experimentally- 
Determined  Apparent  Resistivity 
Curves . 23 

Price's  Theory  of  Magneto-Telluric 

Methods . 26 

General  Theory.  .  . . 26 

Two-Layer  Case  with  Extremely 

Resistive  Substratum  .  29 

Chapter  III:  Measurement  of  Magneto-Telluric  Signals 

General . 32 

Field  Instrumentation . 33 

Absolute  Coil  Calibration . IfO 


Continued 


Page 

Playback  Circuit . 43 

Records  of  Magneto-Telluric  Signals  .  .  .  43 

Chapter  IV:  Analysis  of  Magneto-Telluric  Data 

Introduction . 5l 

Outline  of  Mathematical  Steps  .  52 

Power  Spectrum  Graphs . 57 

Chapter  V:  Results  and  Interpretation  . 

Calculation  of  Experimental 

Resistivity . 64 

Experimental  Apparent  Resistivity 

Curves . 64 

Five-Layer  Case . 65 

Layer  of  Variable  Conductivity . 68 

Conclusion. . 75 

Appendix:  Section  I . .  77 

Bibliography . 78 

Table  1:  Depth  of  Penetration  Given  in  KM  ......  7 

Table  II:  Electrical  Properties  of  Earth  Materials.  .  •  7 

Table  III:  Apparent  Resistivity  of  the  Sedimentary 

Formations  Obtained  from  Electrical  Logs.  69 

Figure  1:  Location  of  Magneto-Telluric  Field  Stations 

and  Other  Regions  of  Interest .  3 

Figure  2a:  Origin  of  Micropulsations .  10 

Figure  2b:  Location  of  Micropulsations  on  This 

Representation  of  Part  of  the  Electro- 
Magnetic  Spectrum .  10 

Figure  3:  Western  Canada  Sedimentary  Basin  ••••••  12 


Continued  Page 

Figure  l\.:  Magneto- Telluric  Regions  of  Interest . 15 

P’igure  5:  Magneto-Telluric  Standard  Curves . 19 

Figure  6:  Instruments  for  Measuring  Electric  Field.  •  •  34 

Figure  7:  Block  Diagram  and  Photographs  of  Instruments 

for  Measuring  Magnetic  Field.  .  .  37 

Figure  8:  Equivalent  Circuit  for  Magnetic  Field 

Measurement . 38 

Figure  9:  Circuits  for  Calibrating  Magnetic  and 

Electric  Components . •  •  •  43- 

Figure  10a:  Calibration  Curves  for  Channels  5  &  6.  .  .  •  1+2 

Figure  10b:  Calibration  Curves  for  Electric  and  Magnetic 

Components . 1+2 

Figure  11:  Block  Diagram  and  Photograph  of  Playback 

Circuit.  . . 44 

Figure  12:  Bon  Accord  Magneto-Telluric  Records  .....  46 

Figure  13:  Bon  Accord  Magne to-Telluric  Records#  ....  47 

Figure  l4 :  Cooking  Lake  Magneto-Telluric  Records  ....  46 

Figure  15:  Cooking  Lake  Magneto-Telluric  Records  ....  49 

Figure  l6:  Irricana  Magneto-Telluric  Records . 50 

Figure  17a:  Bon  Accord  Power  Spectra  Graphs . 59 

Figure  17b:  Bon  Accord  Power  Spectra  Graphs . 59 

Figure  18a:  Cooking  Lake  Power  Spectra  Graphs . 60 

Figure  l8b:  Cooking  Lake  Power  Spectra  Graphs . 60 

Figure  19a:  Onoway  Power  Spectra  Graphs . 6l 

Figure  19b:  Onoway  Power  Spectra  Graphs . 6l 

Figure  20a:  Kavanagh  Power  Spectra  Graphs . 62 

Figure  20b:  Kavanagh  Power  Spectra  Graphs.  .......  6 2 

Figure  21a:  Irricana  Power  Spectra  Graphs . 63 

Figure  21b:  Irricana  Power  Spectra  Graphs . 63 


Continued 


Page 


Figure  22:  Theoretical  Apparent  Resistivity  Curves 

and  Actual  Earth  Models  -  Leduc  and 
Willingdon . 70 

Figure  23:  Experimental  and  Theoretical  Apparent 

Resistivity  Curves  -  Cooking  Lake  .  •  •  71 

Figure  2l\.:  Experimental  and  Theoretical  Apparent 

Resistivity  Curves  -  Onoway .  72 

Figure  25:  Experimental  and  Theoretical  Apparent 

Resistivity  Curves  -  Bon  Accord  ....  73 

Figure  26:  Experimental  Curves  -  Kavanagh .  714- 

Figure  27:  Experimental  Curves  -  Irricana .  7U- 


1 


CHAPTER  I 

INTRODUCTION 

9 

An  investigation  of  the  subsurface  was  carried  out  In 
central  Alberta  by  using  a  geophysical  exploration  method 
termed  "the  magneto-telluric  method"  (Cagniard,  1953)*  The 
results  obtained  were  compared  with  those  from  geological 
and  other  geophysical  investigations, 

A  brief  description  of  pertinent  subsurface  properties 
will  now  be  given.  Extensive  knowledge  of  the  structure  and 
composition  of  the  top  layer,  which  is  referred  to  as  the 
western  Canada  sedimentary  basin,  has  been  obtained  from 
seismic  records  and  from  well  samples.  Moreover,  the 
resistivities  of  the  sedimentary  formations  In  central 
Alberta  are  obtainable  from  numerous  electrical  logs.  Hence 
a  comparison  of  the  average  resistivity  and  thickness  of 
this  layer  as  obtained  by  the  magneto-telluric  method,  with 
those  from  in  situ  measurements  is  possible.  Underlying  the 
sedimentary  basin,  which  is  2.7  km  thick  in  the  Edmonton 
area,  is  the  Precambrian  Shield.  A  combined  geophysical 
and  petrological  study  of  the  Precambrian  in  central  Alberta 
by  Garland  and  Burwash  (1959)  shows  the  lithology  of  the 
basement  to  vary  over  this  area.  Investigation  of  the 
subsurface  at  great  depths  has  been  carried  out  in  Alberta 
at  two  different  locations  and  by  different  geophysical 
methods.  At  Meanook,  Niblett  and  Sayn-Wittgenstein  (i960). 
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by  employing  the  magne to-telluric  method  have  detected  a 
discontinuity  in  resistivity  at  a  depth  of  80  km.  Seismic 
records  at  Irricana  (Richards  and  Walker,  1959)  indicate 
the  depth  to  the  Mohorovicic  discontinuity  to  be  l\.3>  km. 

Beneath  the  Precambrian  lies  the  mantle,  for  which 
our  knowledge  of  the  composition  and  structure  is  limited 
mainly  to  information  obtainable  from  seismic  records. 

Magneto- telluric  recordings  were  taken  at  Irricana, 

Bon  Accord,  Onoway,  Cooking  Lake  and  Kavanagh  in  the  summer 
and  fall  of  1961.  In  Figure  1  are  shown  the  locations  of 
the  above-mentioned  places  together  with  the  locations  of 
wells  in  Alberta  v/hich  have  penetrated  into  the  basement. 

The  results  obtained  by  the  magneto-telluric  method 
were  found  to  be  consistent  with  those  obtained  by  other 
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A  DEPTH  TOELECTRICAL  DISCONTINUITY  DETERMINED  HERE 

Pig.  1.  Location  of  magneto-telluric  field  stations  and 
other  regions  of  interest. 


CHAPTER  II 


THEORY 


The  Magneto- Telluric  Method  (Cagniard,  1933) 

Outline  of  Method 

Consider  a  region  where  the  layers  are  homogeneous 
and  horizontally  stratified  and  where  plane  elec tromagne tic 
waves  strike  the  surface.  Then,  providing  adjacent  layers 
have  an  appreciable  electrical  resistivity  contrast  and 
providing  the  signal  has  relatively  strong  harmonic 
components  in  the  low  frequency  range  (eg:  0.002  to  10  cps 
range),  then  the  thickness  and  resistivity  of  each  layer  is 
obtainable  by  performing  the  following  operations: 

■  (1)  measure  simultaneously  mutually  perpendicular 
components  (Ex  and  Hy  say)  of  the  electric  and 
magnetic  fields  tangential  to  the  earth’s  surface. 

(2)  process  this  data  to  determine  those  harmonic 
components  which  make  a  substantial  contribution. 

(3)  compute  the  "apparent  resistivity",  pa  ,  for  a 
number  of  values  of  the  period,  T,  in  the  period 
range  of  interest  by  means  of  the  following 
formula : 


.  .  .II  -  1 


The  ratio  FJX/HV  is  called  the  "Intrinsic  impedance 
of  the  earth  for  plane  waves  ( Schelkunoff ,  19^*3)  • 
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The  meaning  of  apparent  resistivity  and  intrinsic 
impedance  will  be  dealt  with  in  a  subsequent 
section* 

(4)  plot  on  transparent  log-log  paper  apparent 
resistivity  values  against  period* 

(5)  compare  this  curve  with  "master  curves"  computed 

from  theoretical  models  and  hence  determine  the 
"best-fit"  curve*  v 

(6)  by  means  of  a  few  simple  operations  (to  be 
described  later) ,  determine  the  thickness  and 
resistivity  of  each  layer* 

The  advantage  of  Cagniard’s  method  over  electrical 
prospecting  methods  are  two-fold.  The  depths  of  investigation 
are  orders  of  magnitude  greater  than  those  dealt  with  in 
the  latter  method;  neither  power  supplies  nor  base  stations 
are  required. 

In  theoretical  work,  many  writers  on  this  subject  have 
employed  electromagnetic  units  and  in  practical  applications, 
"practical"  units.  Hence  this  procedure  will  be  followed. 

A  table  for  conversion  from  em  cgs  to  practical  units  is 
given  in  the  appendix  (section  I). 

Frequency  Range  of  Interest 

This  range  is  from  0*002  to  1  cps.  In  Figure  2(b) 
is  shown  the  location  of  magneto-telluric  signals  in  the 
electromagnetic  spectrum.  The  reason  for  recording  waves 
in  this  particular  frequency  range  can  be  explained  by 
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examining  the  concept  of  "skin  depth”  or  "depth  of 
penetration".  It  designates  the  depth  to  which  an  electro¬ 
magnetic  wave  penetrates  into  a  uniform  conductor  before 
its  amplitude  is  reduced  to  l/e  of  its  surface  value.  The 
skin  depth,  p,  is  given  by  the  familiar  expression  (in  era  cgs 
units ) 


.II  -  2 


T  =  period  in  seconds 
p  s  resistivity  in  era  units 
p  s  penetration  in  eras 

In  Table  I  is  shown  the  skin  depth  for  various  periods 
and  resistivities;  in  Table  II  is  shown  for  typical  materials, 
the  skin  depths  at  two  periods. 

In  this  investigation,  the  region  of  interest  extends 
from  the  surface  (sedimentary)  layers  to  the  top  region  of 
the  mantle.  In  order  to  investigate  the  sedimentary  basin, 
magneto-telluric  signals  v/ith  periods  as  small  as  0.1  seconds 
are  required;  for  mantle  studies,  periods  of  hundreds  of 
seconds  are  required. 


Uniform  Fields 

Investigation  of  telluric  current  records  in  France 
and  in  western  Canada  indicate  these  currents  to  be  uniform 
over  large  areas. 

In  France,  Blavier,  upon  examination  of  the  telluric 

current  records  taken  from  parallel  lines  separated  by  as 

much  as  270  kms,  found  them  to  be  strikingly  similar. 
Schlumberger  and  Kunetz  (19^8)  found  a  correlation  between 
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TABLE  I 


DEPTH  OP  PENETRATION  GIVEN  IN  KM  (After  Cagniard) 


0/T  -  1  sec 

4 

ohm  meters 

3  sec 

10  sec 

30  sec  1  min 

2  min 

0.2 

0.225 

0.390 

0.712 

1.23 

1.71)- 

247 

1 

0.503 

0.872 

1.59 

2.76 

3.90 

5.5l 

5 

1.13 

1.95 

3.56 

6.16 

8.7  2 

12.3 

10 

1.59 

2.76 

5.03 

8.72 

1 2.3 

174 

5o 

3.56 

6.l6 

11.3 

19.5 

27.6 

39.0 

250 

7.95 

13.8 

25.2 

lp.6 

6l.6 

07.2 

1,000 

3.5.9 

27.6 

50.3 

87.2 

123 

174 

5,000 

35.6 

6l.6 

113 

195 

276 

390 

TABLE  II 

ELECTRICAL  PROPERTIES  OP  EARTH  MATERIALS  (After 

X 

Garland) 

Typical 

resistivity. 

T  =  1  , 

Depth  of  penetration"' 

sec  T  r  1  min  T  =  1  day 

Material 

ohra  meters 

km 

km 

km 

Sea-water 

.20 

0.22 

1.7 

4.5 

Unconsolidated  soils 

50.00 

3.55 

27.5 

1040 

Sedimentary  rocks 

10Q0Q 

5.03 

3  8.9 

lkoO 

Igneous  rock 

1  x  1# 

50.3 

389 

1I18OO 

Mantle,  1000  lan  depth 

100 

0.50 

3.9 

147 

‘"‘The  depth  of  penetration  is  the  depth  within  the  material  at  which 
an  electric  or  magnetic  field  of  the  period  shown  would  be  reduced 
to  l/e  of  its  surface  value. 
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telluric  records  taken  in  locations  as  far  separated  as 
France  and  Madagascar.  Garland  and  Webster  (i960)  infer 
from  telluric  records  in  western  Canada  that  source  distance 
effects  are  relatively  unimportant.  Simultaneous  records  of 
the  E-IV  component  of  the  electric  signal  taken  in  1961  at 
Edmonton  and  Clive  (N-S  separation  distance  75  miles)  by  the 
Universities  of  Alberta,  British  Columbia  and  California 
show  a  remarkable  correlation. 

However,  at  Meanook,  investigation  of  magneto- telluric 
records  by  Niblett  and  Sayn-Wittgens tein  (1961)  indicate 
that  source  effects  may  be  important  for  that  location. 

For  this  investigation,  interpretation  of  the  results 
will  be  made  by  assuming  that  the  waves  are  uniform  over 
lateral  distances  of  at  least  50  kilometers  and  that  the 
material  in  each  layer  is  isotropic. 

Then  the  electric  (E)  and  the  magnetic  (H)  components 
are  orthogonal  and  consequently  obey  the  following: 


x  =  N-S  component 
y  =  E-’.V  component 


However,  if  the  results  indicate  that  the  above  relations 
do  not  apply,  then  source  and  anisotropic  effects  will  be 
considered . 

Classif ic ation  and  Or igin  of  hicropulsations 

Benioff  (i960)  has  observed  that  magneto- telluric  records 
usually  occur  in  a  variety  of  distinct  types,  each  with  its 

In  Figure  2b  are  shown  pictorial  representations 


own  features. 
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of  the  four  main  types  (A,  B,  C,  &  D) . 

Type  A  Waves 

These  have  periods  up  to  several  seconds  and  often 
persist  up  to  1|8  hours.  The  amplitudes  of  these  waves  are 
often  modulated  by  longer  period  waves  and  are  then  described 
as  "pearly’' . 

Type  B  Waves 

These  have  a  period  from  2  to  10  seconds.  They  are 
connected  v/ith  auroral  activity  and  are  related  to  magnetic 
activity  on  the  earth’s  dark  side. 

Type  C  Waves 

These  range  in  period  from  8  to  lj.0  seconds.  Since  they 
appear  only  on  the  sunlit  side,  they  are  called  ,’day-time,, 
activity.  It  is  thought  that  this  type  modulates  Type  A  to 
form  pearls. 

Type  D  Waves 

The  period  range  is  from  40  to  several  hundred  seconds. 
On  occasion,  long  trains  of  these  waves  have  been  observed; 
at  other  times,  they  occur  as  a  single  pulse. 

In  Figure  2(a)  are  shown  two  different  ways  in  which 
micropulsations  are  thought  to  originate. 

Other  means  have  been  observed  by  which  low  frequency 
waves  are  generated.  For  example,  waves  with  frequencies 
around  10  cps  are  the  resonant  frequencies  of  the  electro¬ 
magnetic  cavity  formed  between  the  conducting  earth  and  the 
ionosphere.  Lightning  probably  starts  this  cavity  ringing 
(Balser  and  Wagner,  i960). 
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ORIGIN  OF  MICROPl’LSATIONS  may  he  far  oul  in  earth’s  mag¬ 
netic  field  or  near  the  earth  in  the  ionosphere.  At  top.  wave  gener¬ 
ated  far  away  by  protons  from  the  sun  mo\e*  toward  earth  as  a 

Fig. 


hydromagnetir  wave,  traverses  ionosphere  <  color),  change*  to  long 
electromagnetic  wuve  t  *oliti  line)  and  strikes  earth.  \l  bottom, 
solar  protons  come  in  close,  generating  long  wave  in  the  ionosphere. 

2(a) 
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Location  of  micropulsations  on  this  representation  of  part 
of  the  electromagnetic  spectrum  (after  Heirtzler,  1962). 


Fig.  2(b) 
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Hence  the  components  of  micropulsauions  have  diverse 
origins  and  the  aetermina tion  of  their  origins  is,  in  general, 
complicated  by  the  presence  of  the  ionosphere. 

Layering  of  Karth 

Cagniard’s  method  is  especially  suited  to  regions  of 
horizontal  layering  of  fairly  homogeneous,  isotropic  properties. 

In  Figure  3  is  reproduced  a  geological  cross  section  of 
the  western  Canada  sedimentary  basin.  This  cross  section 
commences  at  Nordegg  (see  Figure  1)  which  is  in  the  Alberta 
foothills  and  passes  through  Leauc  (20  miles  south  of  Edmonton)  • 
This  line  is  approximately  perpendicular  lo  the  Rocky  Mountains. 
Taking  into  account  the  fact  that  the  vertical  scale  exagger¬ 
ation  is  approximately  lj.2  times,  the  layering  can  be  seen  to 
be  approximately  horizontal  and  of  fairly  uniform  thickness 
over  large  lateral  distances. 

Examinations  of  electric  logs  taken  at  numerous  locations 
in  central  Alberta  indicate  that  the  layering  in  the 
sedimentary  basin  is,  in  general,  electrically  inhomogeneous. 
That  is,  the  apparent  or  mean  resistivity  for  a  given  layer 
was  found  to  vary  at  different  locations.  However,  over 
regions  of  horizontal  extent  of  say  $0  to  80  miles,  this 
change  is  probably  negligible.  Hence,  for  apparent 
resistivity  calculations,  the  layers  were  assumed  to  be 
homogeneous . 

From  geological  investigations,  the  sedimentary  basin 
is  thought  to  be  horizontally  isotropic.  Anistropic  media 
are  observed  mainly  in  shales,  which  often  have  mineral 


RIMBEY  -MORINVILLE 
Reef  Choin 


s 


s 


13 


crystals  of  fairly  elongated  shape  (mica,  kaolin,  etc) 
embedded  in  them.  An  electric  current  travels  with  greater 
facility  along  the  interstices  which  may  exhibit  a  parallel¬ 
ism  with  the  stratification  and  which  are  usually  saturated 
with  the  mineralized  water.  An  experiment  to  test  the  ground 
(to  a  depth  of  100  feet)  for  anisotropy  was  performed  at 
Kavanagh.  The  results  indicate  the  ground  to  be  horizontally 
isotropic  do wn  to  this  depth. 

The  Precambrian  in  central  Alberta  has  been  investigated 
by  Garland  and  Burwash  (i960)  in  a  combined  gravitational 

and  petrological  study.  Their  results  indicate  the  basement 

lithology  to  vary  over  this  area  and  to  consist  of  three  main 
types;  namely,  gneissic,  granitic  and  basic  phases.  Hence 

it  would  be  incorrect  to  assume,  a  priori,  that  the  electrical 
resistivity  of  the  basement  be  the  same  at  trie  various 
recording  sites.  Moreover,  since  relatively  little  is  known 

of  the  electrical  properties  in  situ  of  the  basement,  it  is 
not  known  whether  the  basement  is  isotropic  or  anisotropic. 

Since  micropulsations  penetrate  readily  into  the  basement,  the 

magneto- tellur ic  method  should  be  ideal  for  investigating  the 

electrical  properties  of  this  region. 

Very  little  is  known  about  the  electrical  properties  of  the 

top  portion  of  the  mantle.  Prom  an  analysis  of  rnagne to-telluric 
records  taken  at  Meanook,  Niblett  and  Sayn-Wittgenstein  (1961) 

suggest  that  there  may  be  an  electrical  discontinuity  at  a 
depth  of  80  km. 

Thus  the  subsurface  of  central  Alberta  probably  consists 
of  a  highly  resistive  precambrian  sandwiched  between  two 
relatively  conductive  layers. 
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Displacement  Current  In  a  Conductor 

It  is  evident  from  an  inspection  of  Maxwell’s  relation: 


V  x  M  =  4'"  J  *  D  , 


•  •  •  X  X  —  J 1 8 


that  in  a  conductor  of  conductivity  cr  and  dielectric 
constant  e  ,  there  is  both  a  conduction  current  and  a 
displacement  current.  The  right  hand  side  of  II  -  4a  can 
be  written  (for  sinusoidal  waves)  as: 


.  .  .II  -  Ij-b 


The  displacement  current  term  becomes  less  important  as  the 
wave  frequency  is  decreased.  Figure  Ij.  shows  the  magneto- 
telluric  region  (in  so  far  as  this  investigation  is 
concerned)  and  the  region  where  the  displacement  current 
becomes  appreciable  In  comparison  with  the  conduction  currents. 
Examination  of  the  figure  shows  that  in  the  region  of  interest, 
the  displacement  current  is  negligible.  Henceforth,  the 
displacement  current  term,  when  dealing  with  the  earth,  will 
be  omitted. 

Direction  of  Propagation  of  Magneto-Telluric  Wave  in  the  Earth 

Consider  a  plane  wave  of  propagation  constant  k-^ 
striking  the  earth  at  an  angle  of  incidence  <3(  .  Let  ko  and 
Oz  be  the  corresponding  value  for  the  refracted  ray.  From 
Snell’s  law: 


.  .n-5 


RESISTIVITY  (ohm  meters) 
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10  10  10' 

FREQUENCY  (cps.) 


Pig.  1|..  Location  of  region  where  displacement  current 

becomes  larger  than  0.1$  of  conduction  currents 
(after  Neves*  1957)  and  magneto-telluric  regions 

of  interest. 


l6 


The  wave  equation  is: 


V  H 


(4  ii  i  U  cr  ~  u)2e  )  H  .  ....11  -  6 


2.  2. 

Hence  K,  =  -  coLe  and  Kz  =  4'ftLbj(T  .  .  .  .11-7 

|3,n  f  I -  O 

since  the  conduction  current  term  (on  the  denominator)  is 
much  greater  than  the  displaceraent  current  term  (on  the 
numerator)  in  the  magne to- telluric  region  of  interest.  An 
inspection  of  Figure  lj_  shows  this.  Thus,  no  matter  what 
the  angle  of  incidence,  the  direction  of  propagation  in  the 
earth  will  be  essentially  vertically  downward. 


Intrinsic  Impedance  of  the  Earth  which  is  Treated  as  a 

Semi-Infinite,  Uniform  Conductor 

Cartesian  co-ordinates  with  the  X-Y  plane  on  the  earth* s 
surface  and  z  pointing  downwards  will  be  used.  X  points 
towards  the  geographic  north  pole. 

Consider  a  plane  wave  with  the  electric  vector  polarized 
in  the  N-S  direction  and  the  magnetic  vector,  in  the  E-W 
direction,  impinging  upon  the  earth.  From  Maxwell’s  curl 
relation  II  -  Ij.,  one  obtains  (for  the  conductor): 


Ihr  =  -  4-  7f  cr  £  x 

3  *  •  •  • 

From  the  wave  equation  II  -  6,  one  obtains  (for  the 

conductor) : 


II 


/4  n  k  cr 


8 


y 


.II  -  9 


17 


The  Intrinsic  impedance  of  a  uniform  earth  for  plane  waves  is 


defined  as  the  ratio  of  Ex  to  Hy  where: 


(ZcrTJ 


"  Vz 


.  .11-10 


Intrinsic  impedance  is  analogous  to  the  "characteristic 
impedance"  of  a  one-dimensional  transmission  line.  For 
free  space,  this  impedance  is  376.6  ohms. 

The  resistivity  p  can  be  expressed  in  terms  of  the 
intrinsic  impedance  and  is  given  by  II  -  1. 

After  a  brief  discussion  on  apparent  resistivity,  it 
will  be  shown  how  equation  II  -  1  is  used  to  interpret 
subsurface  structure. 

Apparent  Resistivity 

Consider  two  models:  a  Model  (A)  which  is  horizontally 
stratified  into  n  electrically  homogeneous  layers  of 
resistivities  f  p2  t  • • •  pn  >  and  also  a  Model  (B)  of 
uniform  resistivity  p A  •  Consider  a  wave  of  period  T  which 
penetrates  into  the  nrn  layer  of  Model  (A);  let  the 
intrinsic  impedance  for  this  Model  be  Ey/ily  l%_o  •  Not/ 
let  a  wave  with  the  same  period  T  impinge  upon  Model  B) . 

Let  be  of  such  a  magnitude  that  the  intrinsic 
impedance  for  Model  (B)  is  the  same  as  that  for  Model  (A). 
Then  pA  is  called  "apparent  resistivity"  for  the  n  layers 
(for  that  particular  value  of  T)  and,  physically,  it  is  a 

measure  of  the  mean  or  average  of  the  resistivities  of  the 

n  layers. 
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The  Two-layer  Earth  Model 

At  certain  locations,  (eg:  where  the  Precambrian  Shield 
extends  from  the  surface  to  the  mantle)  the  earth  approximates 
a  two-layer  model.  In  Figure  5(a)  are  shown  graphs  of  apparent 
resistivity  versus  period  (on  log-log  paper)  for  a  number  of 
resistivity  values  of  the  bottom  layer,  the  upper  layer 
being  1  kilometer  thick  and  having  a  resistivity  of  1  ohm 
meter. 

The  expression  for  the  intrinsic  impedance  for  plane 
waves,  Ex/l-Iy  j  4  o  is  readily  obtainable  for  the  two-layer 
case  from  the  five-layer  case  given  by  II  -  18  by  putting 

P5  =  Pij.  =  P3  =  P 2  and  hjp  5  h3  ■  h2m  hx. 

Two  special  case  of  interest  will  now  be  treated: 

( a )  the  extremely  resistive  substratum  case . 


becomes : 


v  I 


cosK  (*yPj  +  cos  (*yp, ) 

cosK  OVp.)  -  cos  (»h/p.)  ,  .  .  .  .II  -  11 


where  p^  Is  the  skin  depth  in  the  first  medium 
and  p is  its  resistivity. 

When  h  <<  p^,  II  -  11  becomes: 


.  .  .11-12 


The  corresponding  expression  for  p d  is: 


.11-13 


P  A  fL  M 


(a ) 


Two-layer  standard 


curve  s . 


I  _0  10  100  1000  10,000 

T  (SEC) 


(b)  Three-layer  standard  curves. 


Fig.  5.  Magneto-telluric  standard  curves  (after  Yungdl,  19^1)- 


An  interesting  feature  of  the  graph  of  ^  versus 
T  (when  plotted  on  log-log  paper)  is  observed  for 
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this  case.  As  T  increases,  the  graph  of  versus 
T  becomes  linear  and  its  slope  makes  an  angle  of 
Jpf?0  with  the  abscissa.  To  see  this,  take  the  log 
of  both  sides  of  II  -  13.  One  gets: 

Joy  pa  =  Jo<f  t  77?W  +  ^7  1  ;  .  .  .xi  .  ii,. 

Taking  the  derivative  of  both  sides  with  respect 


to  log  T,  one  gets: 

d 


=  / 


ci  ( Joy  T)  ^  //d/  1  •  ....II 

(b)  the  extremely  conductive  substratum  case. 

when  P2  =  0,  the  expression  for  E^/Hy  becomes: 


-  15  \ 


The  corresponding  graph  of  pA  versus  T  for  large 
values  of  T  becomes  linear  and  has  a  slope  of  -if5°. 

For  the  cases  where  the  experimental  apparent 
resistivity  curves  (when  plotted  on  log-log  paper) 
have  slopes  with  magnitudes  greater  than  i|-50  with 
respect  to  the  abscissa,  interpretation  of  sub¬ 
surface  structure  by  Cagniard’s  method  is  not 
possible.  Then  such  factors  as  source  effects  and 
anisotropy  must  be  taken  into  account.  However,  it 
is  only  when  these  two  effects  can  be  dealt  with 
individually  that  an  interpretation  becomes  possible. 
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Five -Layer  Model 

An  Inspection  of  Figure  3  shows  the  sedimentary  basin  in 
the  Leduc  area  to  consist  of  rocks  from  five  geological  time 
divisions;  they  are  the  UPPER  CRETACEOUS,  the  LOWER  CRETACEOUS, 
the  UPPER  DEVONIAN,  the  MIDDLE  DEVONIAN  and  the  UPPER  and 
MIDDLE  CAMBRIAN  periods.  In  order  to  determine  the  influence 
these  five  layers  have  on  the  shape  of  the  apparent  resistivity 
curve,  expressions  for  the  intrinsic  impedance  and  the  corr¬ 
esponding  apparent  resistivity  were  developed  for  a  five- 
layer  case.  The  following  expressions  (which  are  original 
calculations)  were  derived  because  previous  calculations  , 
were  for  two  or  three  layers  only. 


|’T  LAYER 

2"‘°LA  YE  R 

3*aLAYER 

p 

4™  LAYER 

5rM  LAYER 
i 

FIVE  -  LAYER 
EARTH  SECTION 


Then : 

P« 


Let: 

k  -  subscript  denoting  layer  no, 
h^  =  depth  to  bottom  of  layer  k 
p ^  =  resistivity  of  layer  k 

T  *  period  of  wave 
pi^  -  skin  depth  relative  to  layer  k 

=  . 

p  z  apparent  resistivity 


II  -  18 


where  : 


E 

E 


R 


I 


H 


R 


H 


I 


D2j.,c  o  s  h  QtJ_, 
D2j  ,s  i  n  h  9ZJ_, 
D;„s  i  n  h  9lj., 


Ju  L 


t  U.,cosh  dtJ., 


+ 

+ 

+ 

+ 


(-1  /”D„  sinh 6L 
(-1  f  D2J  cosh  6ZJ_, 
(-1)  D,cosh  9ZJ_, 
(-1  )  U  sinh  Qu_t 


•COS 9U_, 

•sin  etJ_, 

■  COS  9ZJ., 
sin  9U., 


i 


(cont'd  J 
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and 


D, 

D. 

D 

D. 

D. 

D. 


D 

D 

D 

D, 

D 

f 

D 

D 

D 
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(i+jrX^X^X^) 

(?+W 

I  2  2  3  3  4  4 

(i+ft^WX*) 

12  2  3  3  4  r5 


I  2  2  3  3  4  4 


D7  -  (£+pXp-pXp+fXp) 

'  rl  ~2  ~2  r3  r3  ^  r4 


(i+iX^-^X-p+^(^) 

1  I  2  2  '  3  4  5 


(tmmv 

(X-pXp-XXjj+^(  p) 
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(i-p)(-p-i)(-p-?)(i) 

I  2  2  3  3  4  4 

/  _L  _  iXLi  _  IVi  __  J_V_L  ) 

'  P  p/^P  p/'P  PA  P  ' 

r,  f~2  ^2  r3  ^  4  r5 


(  jj  _  ^Xp+  p )(p"  (5X  p5 ) 


( p  pXp+pXp+pX  p ) 


$,  -  h,(^-q)+h(^-^}+^(p3-^)+h,(^) 


$3  =  MrW-pH^-^1 

rl  2  2  r3  K3  *4  *4 


$5  _  hX^-pl+h^  +  ^-h^+^+hXp  ) 

I  r2  2  2  *3  3  ^3  4  r4 


@7  -  h,(p,  py\(p2+ p)  hjp3  p4)  h(p 


$9  h,(p+J  hJp  +  p)+hip,  pJ+hX  pj 


I  T  2 


2  '3  “ ' 3  '4 


$n  =  h(p  +  p)-h2(p  +  p)+h3(p+p)-h4(p  ) 

12  2  3  3  4  4 


$13  hX  pf  +  p2)  hX  p,  p  X  hi  P3  +  p}  p4  ) 


$i5  =  h,  (p  +  pWp  -  pHtp-  p)~  hi -p ) 

I  r2  2  3  3  ^  ^ 
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Interpretation  of  Experimentally-Determined  Apparent 

ResisTlvity  Curves" 

Interpretation  of  experimental  resistivity  curves  in  terms 
of  layering  In  the  subsurface  can  be  effected  by  employing 
"standard11  curves  such  as  are  shown  in  Figure  Since  this 

method  requires  the  use  of  the  law  of  elec trodynamic  simili¬ 
tude,  a  brief  derivation  of  this  law  will  be  given. 

Suppose  the  two  models  shown  below  are  connected  by  the 
following  relations: 


hr  ^ 

ft. 

Earth  Model 

ft. 

^IS  fis 

ft. 

J 

P 2.S 

71 

Standard  Model 

=  K'/l  -  i 

--  ^pti  -  n 

=  K Lhls  -  hi 

-  KtT5  -  iv 

.  .  .  .11-19 


p  -  resistivity 
h  ■  thickness 
T  =  period 

K  s  dimensionless  measure  number 


Consider  two  waves  with  equal  amplitudes  impinging  upon  both 
models.  Then  the  law  of  similitude  stipulates  that  the 
wave  equation  be  satisfied  for  both  models.  That  is: 


+ 


(a) 

.11  -20 


0 
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and  also: 


-  /4 n(^K^)E 


0 


(b) 

.  .  .  .II  -  20 


By  using  equation  (III),  II  -  19,  one  obtains: 


.11-21 


Hence  by  using  relations  in  II  -  19,  equation  (a)  of 
II  -  20  can  be  transformed  into  the  following  form: 


-fev,  E 


4Tr^,X^jE 


-  0 


or 


V/E  <  - 


=  0 


.  .  .  .II  -  22 


This  relation  becomes  identical  to  equation  (b)  of 
11-20  when: 

K L  =  Kt  KP  _  ...  .ii  -  23 

This  then  is  the  condition  which  the  dimensionless  measure 
numbers  must  satisfy  if  the  lav/  of  similitude  is  to  hold. 

Suppose  now  that  one  has  an  experimental  apparent 
resistivity  curve  for  a  two-layered  earth  and  one  wishes  to 
determine  its  parameters.  Let  the  parameters  be  represented 
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by  those  of  the  earth  model;  let  the  standard  model  have 
p-^  =  1  unit  and  h-^  r  1  unit.  Then  the  procedure  for  deter¬ 
mining  these  parameters  is  as  follows: 

(1)  place  this  experimental  apparent  resistivity  curve 
(drawn  on  transparent  paper)  over  the  standard 
curves  (Figure  5a)  and  translate  it  until  a  good 
fit  is  obtained. 

(2)  to  obtain  p^9  read  the  value  of  the  ordinate  on 
the  top  graph  which  coincides  with  the  curve 

p  z  1  on  the  bottom  graph. 

(3)  to  obtain  t  read  the  value  of  pxS  which  best  fits 
the  experimental  curve  and  multiply  this  value 
with  n  . 

I  J  XL 

(Ij.)  to  determine  hlA,  read  the  value  of  the  abscissa 
(  T,4)  on  the  top  graph  which  coincides  with  TJS  s  1 
on  the  bottom  graph.  Then  from  equations  (III)  of 
II  -  19  and  II  -  23,  one  gets: 

h  ljt  —  Kuhts  —  V  Kp  h ti 

Since  h  =  l  and  plS  =  l  , 


h\  x  ~~  and  Kp  Pijl. 


.  .  .  .II  -  2if 


Hence  an  Interpretation  is  obtained. 

When  the  shape  of  the  experimental  apparent  resistivity 
curve  was  suggestive  of  a  two-layer  case,  the  above  method 
was  used  for  an  interpretation.  However,  many  of  the  curves 
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were  suggestive  of  three  and  four-layer  cases;  for  such 
instances,  theoretical  apparent  resistivity  curves  were 
obtained  by  using  the  formula  for  the  five-layer  case 
(II  -  10) .  A  computer  program  was  written  to  determine 
theoretical  apparent  resistivity  curves  for  layerings  up 
to  five  in  number  and  for  any  resistivity  combinations. 

Pricers  Theory  of  J.Iagneto  -telluric  Methods  .Then  the  Source 
is  Considered 

Although  Price1 s  theory  (1962)  is  not  used  to  obtain  a 
quantitative  interpretation  of  the  experimental  apparent 
resistivity  curves  obtained  in  this  investigation,  it  is, 
nevertheless,  suggested  as  a  possible  means  of  interpreting 
features  of  these  curves  that  are  difficult  to  explain  by 
Cagniard’s  method.  Consequently  a  brief  description  of 
Price fs  theory  will  be  given. 

Price  took  into  account  the  nature  of  the  source 
field  in  a  general  way.  His  method  includes  Cagniard’s 
plane-wave  treatment  and  Wait*s  (195^1-)  dipolar  source  as 
special  cases. 

General  Theory 

Cartesian  co-ordinates  oriented  with  the  X- v  plane  on 
the  earth1 s  surface  and  Z  vertically  downwards  will  be 
chosen.  For  the  present  the  earth  will  be  treated  as  a 
semi-infinite  conductor. 

An  expression  for  the  intrinsic  impedance  ^x/ry/e.,0 
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of  the  earth  when  the  lateral  dimensions  of  the  source  field 


are  considered  will  be  derived  starting  from  Maxwell’s 
equations.  Neglecting  displacement  currents  in  the  conductor 
and  taking  unity  for  the  permeability  of  the  earth.  Maxwell’s 
equation  becomes: 


V  X  tt  ~  477  l  =■  4- IT  c r($)E 


.  .II  -  25 


V  x  E  =  -  i  u  H  . 


.  .II  -  2 6 


The  wave  equation  obtained  from  the  above  equations  in  the 
usual  manner  is: 


.  .II  -  27 


In  a  previous  paper  (Lahiri  and  Price  1939) t  it  was 
shown  that  currents  induced  in  a  conductor  (where  layers 
of  equal  conductivity  are  parallel  to  the  surface)  by  an 
external  field  flow  parallel  to  the  surface.  Although  there 
is  disagreement  on  this  point,  if  one  does  accept  it  as  being 
correct,  then  one  obtains  solutions  to  the  wave  equation  of 
the  form: 


.  .11-28 


where : 


.  .11-29 


Substituting  II  -  20  and  II  -  29  into  the  wave  equation 
II  -  27,  and  employing  the  method  of  Separation  of 
Variables,  one  gets: 
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3  x2  3  $ 


O 

j 


II  -  30 


and 


II  -  31 


Let: 


0 


2. 


2. 


+  4 TT  i  uj  <T(?  ) 


II  -  32 


^  is  a  constant  and  the  reciprocal  of  u  is  shown  by 
Price  to  be  a  measure  of  the  horizontal  extent  of  the  source 
field.  The  range  of  the  values  for  is  may  be  obtained  by 
considering  ir/is  as  a  measure  of  the  linear  dimension  of  the 
source  field.  The  least  value  for  v  is  obtained  when  the 
source  dimension  is  equal  to  the  circumference  of  the  earth. 
This  case  corresponds  to  the  longest  wavelength  of  the  field. 
To  obtain  the  greatest  value  for  V  of  practical  importance. 
Price  considers  the  case  where  the  linear  dimensions  of  the 
source  are  four  times  the  height  (100  kms)  of  the  ionospheric 
jet  currents.  The  range  of  v  thus  obtained  is  1*5?  x  10  ^ 
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to  1.57  X  10  ^  cm"*^* •  Cagniard’s  case  corresponds  to  taking 


v  equal  to  zero.  Physically,  this  means  that  the  lateral 
dimension  of  the  source  field  is  infinite;  hence  one  is 
dealing  with  plane  waves  at  the  earth’s  surface. 

The  magnetic  field  obtained  by  using  II  -  28,  II  -  29, 
and  II  -  26  is: 


.  .11-33 


The  complex  ratio  Ex/lly  is  obtained  by  using 
II  -  28,  II  -  29  and  II  -  33,  and  is: 


...  .II  -  34- 

Two-layer  Case  with  Extremely  Resistive  Substratum 


By  employing  relation  II  -  3^4- 
and  the  boundary  conditions 


(Ti-p 


that  the  tangential  components 
of  E  and  H  are  continuous  at 


cr  =  0 


z  3  0  and  h,  one  obtains  for  the 
intrinsic  impedance  the 
following: 


11  -  35 


where : 


2 


6 


-  y  +  4  77  i  i-JcT 


By  putting  i>  s  0  in  II  -  35>  &  result  equivalent  to 
Gagniard ' s  II  -  11  is  obtained.  In  this  case,  one  is 
dealing  with  plane  waves. 

'-hen  h-^oo  ,  11-35  becomes: 
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when  ^  is  small,  II  -  36  becomes: 
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.II  -  37 


Waitfs  expression  for  this  case  is: 


Ex)  -  /HI  +  ,  ?  $10* )  ,  o_  ( y-* 

Hy  L  "  v  <*•»•«■  2r*Ws  y  3  ^  3  x1 

...  .II  -  38 

where  yi=4'uci o<r  is  the  square  of  the  intrinsic  propagation 
constant  of  the  earth.  Wait  obtains  the  above  expression 
by  considering  the  field  at  the  earth’s  surface  to  be 
composed  of  a  superposition  of  elementary  waves.  His 
expression  II  -  38  reduces  to  that  of  Cagniard’s  II  -  10 
when  the  second  order  space  derivative  of  the  components 
of  the  magnetic  field  are  small  enough  so  that  the  second 
term  on  the  right  of  II  -  38  can  be  neglected  in  comparison 


31 


with  the  first.  Wait’s  expression  can  be  seen  to  correspond 
to  II  -  37  by  using  II  -  33,  II  -  30  and  II  -  31.  The 
ratio  Kx/ (j  Hy  varies  as  V  is  increased  from  0  to  10“^.  The 
larger  V  is,  the  smaller  Ex/jHy  is  and  also  the  apparent 
resistivity 

An  inspection  of  II  -  32  shows  that  the  relative  importance 
of  u  is  greater  the  smaller  the  values  of  uj  on d  cr  .  On  the 
prairies,  where  the  sedimentary  strata  have  fairly  high 
conductivity  values,  departure  of  the  pd  vs  T  curves  from 
those  obtained  by  employing  Cagniard’s  theory  should  be 
negligible  for  small  periods.  However  at  large  periods,  the 
departure  may  be  significant.  In  places  where  the  surface 
layer  has  a  very  small  conductivity,  such  as  the  area  to  the 
east  of  the  prairies,  departure  of  the  curves  from 

Cagniard’s  should  occur  at  shorter  periods* 
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CHAPTER  III 

MEASUREMENT  OF  MAGNETO- TELLURIC  SIGNALS 


General 

Measurements  were  made  on  the  earth’s  surface  of 
horizontal  components  of  the  total  electric  and  magnetic 
fields.  The  geographic  north-south  and  east-west  components 
of  these  fields  were  measured.  A  component  of  the  electrical 
potential  gradient  was  determined  by  first  measuring  the 
voltage  between  two  electrodes  (driven  into  the  ground)  and 
then  by  dividing  this  voltage  by  the  electrode  separation 
distance  (1,000  feet  in  this  case).  The  corresponding 
magnetic  component  was  determined  by  measuring  the  voltage 
induced  in  a  coil  by  the  varying  magnetic  field. 

Two  special  properties  of  magneto- telluric  signals  are 
their  very  low  frequencies  and  minute  amplitudes.  Because 
of  the  latter  characteristic,  these  signals  are  also  referred 
to  as  micropulsations.  Electric  fields  are  of  the  order  of 
a  hundred  microvolts/kilometer;  magnetic  fields,  of  the  order 
of  a  tenth  of  a  gamma,  where  1  gamma  {X  )  =  10  gauss  (  P  )  . 
The  intensity  of  the  earth’s  field  is  in  the  range  0.25  to 
0.75  gauss;  thus  the  magnetic  signals  of  interest  are 
approximately  10“°  times  that  of  the  earth's  main  field. 

Amplifiers  used  to  record  these  micropulsations  should 
have  the  following  properties: 

(1)  Large  Gain 

Since  the  signals  are  recorded  on  magnetic  tape 
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which  has  a  minimum  input  level  of  1  volt  RMS, 
a  gain  in  the  order  of  thousands  is  required  for 
the  electric  signal  and  hundreds  of  thousands  for 
the  magnetic  signal. 

(2)  Low  noise-level  over  the  frequency  range  of  interest 

Components  of  micropulsations  with  frequencies 
greater  than  0.2  cps  have,  in  general,  extremely 
small  amplitudes.  Hence,  if  an  interpretation 
over  the  frequency  range  0,002  to  1  cps  is  to  be 
possible,  the  noise  level  must  be  very  low  in  the 
0.2  to  1  cps  range. 

(3)  Stability  at  these  long  periods 

The  reason  for  recording  the  signals  on  magnetic  tape 
is  that  the  data  could  be  analyzed  entirely  by  electronic 
means,  thereby  eliminating  what  has  hitherto  been  the  most 
time-consuming  part  of  the  whole  operation. 

Field  Instrumentation 

(1)  Detection  of  electric  signal 

A  block  diagram  of  the  equipment  is  shown  in 
Figure  6. 

The  electrodes  were  of  copper  tubing  2  feet  in 
length  and  lj  inches  in  diameter.  The  electrode 
spacing  was  1,000  feet.  This  type  of  electrode 
was  chosen  because: 

(a)  a  good  electrical  contact  with  the  ground 
could  be  obtained  by  simply  driving  the  tubes  into 
the  ground. 


311- 


electrodes 


DC  REJECTION  FILTER 


LOW  -  PASS  6  db 

TT  F I LT  E  R  AT  2 


AMPLIFIER 


GAIN 


60  cps.  REJECT  FILTER 
TWIN-T  CIRCUIT 


FM  MAGNETIC  TAPE  RECORDER 


REJECTION 

cps. 


xIO  to  x  10,000 


Fig.  6.  Instruments  for  measuring  Electric  Field. 
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(b)  although  the  electromotive  force  between 
copper  and  a  solution  of  its  salt  is  relatively 
large  (-0.60  volts),  changes  in  contact  potential 
are  negligible  because- of  the  short  time  interval 
(I4.O  to  60  minutes)  during  v/hich  recordings  were 
made . 

A  220^f  capacitor  was  inserted  to  block  off  any 
DC  potential  appearing  across  the  electrodes.  On 
the  prairies,  this  potential  is  generally  of  the 
order  of  a  few  millivolts  for  an  electrode  spacing 
of  1,000  feet.  The  low-pass  filter  was  inserted 
to  eliminate  the  high-frequency  components. 

The  amplifier,  which  received  power  from  a  ll5v, 
60  cps  generator,  has  a  maximum  gain  of  10,000. 
However,  above  a  gain  of  1,000,  these  amplifiers 
became  unstable.  Hence  the  gain  was  set  at 
1,000  for  nearly  all  the  records  taken  in  the 
summer  and  fall  of  1961. 

Next  a  twin-T  circuit  with  a  zero  transfer 
admittance  at  60  cps  was  inserted  to  reject  the 
60  cps  chopper  frequency  introduced  by  the 
amplifier.  The  characteristic  impedance  of  the 
filter  was  10,000  ohms. 

Two  devices  were  used  to  monitor  the  signals: 
they  were  a  paper  recorder  (two-channel  BRUSH 
RECORDER  MARK  II)  and  an  oscilloscope  (TECTR0NIX 
TYPE  502  DUAL-BEAM) .  When  the  signal  was  observed 


36 


to  have  long  period  components  of  sufficient 
amplitudes,  it  was  recorded  on  magnetic  tape 
( PRECISION  INSTRUMENT,  PS-200  RECORDER).  This 
frequency-modulated,  7-track  record/reproducer 
has  an  error  of  i  J-db.  The  total  harmonic 
distortion  over  the  frequency  range  of  interest 
is  1.5$.  The  input  signal  range  is  1.0 
volt  RMS;  the  input  impedance,  10,000  ohms. 

(2)  Detection  of  magnetic  signal 

A  block  diagram,  together  with  photographs  of 
some  of  the  instruments  is  shown  in  Figure  7.  The 
equivalent  circuit  is  shown  in  the  following  figure. 

Each  coil  was  made  of  30,000  turns  of  //2 6  wire 
wound  around  a  highly  permeable  core (cone tic  -  A.A.). 
To  lessen  wind  noise,  each  coil  was  put  in  a 
shallow  trench  and  covered  with  soil.  A  30^^ 
capacitor  was  inserted  in  parallel  with  the  coil; 
the  resulting  coil  resonant  frequency  was  1  cps. 

The  voltage  induced  in  each  coil  was  generally 
less  than  a  microvolt. 

Next  a  galvanometer  with  a  taut-band  suspension 
and  a  current  sensitivity  of  3  x  10”^  amps/mm  at 
1  meter  was  inserted.  A  collimated  beam  of  light 
was  directed  onto  the  galvanometer  mirror  from 
where  it  reflected  onto  a  narrow  gap  between  two 
silicon  photoresistor  cells.  The  latter  have  the 
special  property  of  having  an  electrical  resistivity 
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Fig.  7 


Block  diagram  (a)  and  photographs  (b)  of 
instruments  for  measuring  MAGNETIC  field 
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Fig.  8.  Equivalent  Circuit  for  Magnetic  Field  Measurement. 
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which  depends  very  markedly  upon  the  amount  of 
light  falling  upon  them.  The  greater  the  amount 
of  light,  the  lower  is  their  resistance.  Hence, as 
can  be  seen  from  Figure  8,  the  potential  at  point 
P  depends  upon  the  current  flowing  through  the 
galvanometer.  This  combination  of  photoresistor 
circuit  and  galvanometer  constitutes  the  second 
stage  of  the  amplification  process,  the  first 
stage  being  that  resulting  from  the  highly 
permeable  core. 

The  signal  undergoes  another  amplification, 
this  time  by  a  DC  amplifier.  Part  of  this  re¬ 
amplified  signal  is  fed  back  to  the  galvanometer 
to  null  the  current  flowing  through  the  latter. 

By  means  of  this  three-stage  amplification 
process,  micropulsations  of  the  order  of  a  tenth 
of  a  gamma  are  amplified  sufficiently  so  that 
they  can  be  recorded  on  magnetic  tape. 

(3)  Calibration  Circuit 

(a)  Electric  calibration  circuit 

In  Figure  9  is  shown  the  electric 
calibration  circuit.  A  KROHN-HITE  oscillator 
(Model  IpipO A )  was  set  at  1  volt  RMS  and  the 
frequency  was  set  at  certain  values  in  the  range 
0.01  to  10  cps.  Two  recordings  on  magnetic  tape 
were  made  of  this  signal.  One  was  the  signal 


appearing  at  the  output  of  the  oscillator  and  the 
other,  the  signal  after  passing  through  the 
amplifier  circuit.  Because  a  signal  of  1  volt  RMS 
is  much  too  large  for  the  amplifier,  an  attenuator 
was  inserted.  Calibration  curves  (volts  out/milli¬ 
volts  in  versus  frequency)  for  the  electric  circuit 
are  shown  in  Figure  10. 

(b)  Magnetic  calibration  circuit 

This  circuit  is  shown  alongside  the  electri 

calibration  circuit  in  Figure  9.  A  KROHN-HITE 
oscillator  and  a  600  ohm  resistor  were  placed  in 
series  with  the  two-turn  calibration  loop.  A 
calibration  signal  from  this  oscillator  was  sent 
through  the  two-turn  loop  at  a  number  of  frequencies 
in  the  range  0.01  to  10  cps.  The  signal  at  the 
output  of  the  oscillator  was  put  on  magnetic  tape; 
the  calibration  signal  was  also  detected  by  the 
main  coil,  amplified  and  then  recorded  on  magnetic 
tape.  This  calibration  yields  volts  out/milliamp 
in  the  two-turn  loop.  The  calibration  curve  for 
the  magnetic  circuits  is  shown  in  Figure  10b. 

Absolute  Goil  Calibration 

The  above  calibration  circuits  give  volts  out/milliamp 
in.  In  order  to  obtain  volts  out/gamma  in,  (in  practical 
units,  the  magnetic  field  is  measured  in  gammas)  the  box 
containing  the  main  coil  plus  the  two-turn  loop  was  placed  in 
a  solenoid  10  feet  in  length,  1  foot  in  diameter  and  consist¬ 
ing  of  300  turns.  A  signal  from  a  low  frequency  oscillator 
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was  passed  through  the  two-turn  calibration  loop  and  the 
voltage  induced  in  the  solenoid  circuit  was  measured.  The 
magnetic  flux  passing  through  the  solenoid  could  then  be 
computed  and  was  found  to  be  4*5  gamma s/mill lamp  in  the 
two- turn  loop. 

Playback  Circuit 

In  Figure  11  are  shown  a  diagram  of  the  playback  circuit 
together  with  a  picture  of  the  corresponding  apparatus.  The 
data  (in  analog  form)  on  magnetic  tape  were  digitized  and 
recorded  on  paper  tape.  A  Brush  recorder  was  inserted  in 
parallel  with  the  circuit  to  monitor  the  signal  being  punched 
on  paper  tape.  The  SOROBAN  TAPE  PERFORATOR  (Model  GP2) 
records  hole  patterns  at  a  maximum  rate  of  300  characters  per 
second.  The  digitized  data  were  now  in  a  form  whereby  they 
could  be  read  into  the  electronic  computor  (IBM  1620  at  the 
University  of  Alberta)  via  the  paper  tape  reader  (IBM  1621) . 

Records  of  Magneto-Telluric  Signals 

In  Figures  12  to  l6  are  shown  records  taken  in  Central 
Alberta  at  Bon  Accord,  Cooking  Lake  and  Irricana. 

Components  of  magne to-telluric  signals  usually  have  very 
small  amplitudes  above  a  frequency  of  approximately  0.2  cps. 
Most  of  the  records  taken  in  the  summer  and  fall  of  19&1 
showed  a  lack  of  high  frequency  components  in  the  range  0.2  to 
1  cps.  Exceptions  were  those  taken  at  Bon  Accord  (Figures  12 
and  13)  and  Cooking  Lake  (Figure  lip  and  1$)  •  The  magnetic 

components  in  Figure  12  show  clearly  the  presence  of  a  high 
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frequency  component  of  approximately  0.5  cps.  In  Figure  13, 
which  is  a  continuation  of  the  preceding  record,  the  paper 
speed  was  increased  approximately  12  times  to  enhance  the 
presence  of  these  relatively  high  frequency  components. 

The  strongest  signals  recorded  were  at  Cooking  Lake. 

High  frequency  components  can  be  seen  to  be  superimposed 
upon  much  lower  frequency  components  of  considerable 
amplitudes.  Figure  15,  which  is  a  continuation  of  Figure 
lip  but  at  a  paper  speed  10  times  as  high,  shows  more  clearly 
the  presence  of  these  high  frequency  components.  The  records 
taken  at  this  station  yielded  useful  information  over  a  wider 
frequency  range  than  those  taken  at  the  other  four  stations. 

In  Figure  l6  is  shown  a  section  of  the  record  taken 
at  Irricana.  At  this  location,  the  depth  to  the  Mohorovicic 
discontinuity  was  determined  from  seismic  records  (Richard 
and  Walker  1959)*  It  was  hoped  that  magneto- telluric 
measurements  at  this  site  would  yield  information  on 
subsurface  structure  down  to  this  seismic  discontinui ty 
(Jp3  kilometers).  In  order  to  be  able  to  interpret  the  earth's 
structure  to  this  depth  in  central  Alberta,  periods  in  the 
range  of  100  to  300  seconds  are  required.  An  inspection 
of  the  Irricana  record  shows  the  longest  period  present  to 
be  of  the  order  of  100  seconds.  Subsequent  spectral  analysis 
of  this  record  corroborates  the  above  statement. 
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CHAPTER  IV 

ANALYSIS  OF  MAGNETO -TELLURIC  DATA 

Introduction 

Magneto-telluric  signals  are,  in  general,  a  super¬ 
position  of  signals  from  many  different  sources.  Hence 
a  suitable  method  for  determining  the  harmonic  components 
which  make  a  significant  contribution  would  be  a  Fourier 
analysis.  However,  because  of  the  noise  introduced 
through  the  instruments  (amplifiers,  generator,  etc)  and 
external  sources  (power  lines,  lightning,  etc)  a  more 
appropriate  method  would  be  one  in  which  much  of  the  noise 
is  eliminated  (by  statistical  means)  before  the  Fourier 
analysis.  Such  a  method  was  outlined  by  Blackman  and  Tukey, 
(1958)  and  was  used  in  this  thesis. 

Previous  methods  for  analyzing  these  micropulsations 
have  been  the  following: 

(a)  selecting  parts  of  the  record  where  the  signal  is 
very  nearly  sinusoidal  and  measuring  the  amplitude 
and  corresponding  frequency  (Niblett  and  Sayn- 
Wittgenstein,  i960). 

(b)  employing  a  sufficient  number  of  narrow  band-pass 

filters  to  cover  the  frequency  range  of  interest 
(eg:  pass  bands  at  0.005  -  0.02,  0,02  -  0.06, 

etc,  to  0.6  -  1  cps). 

(c)  mathematical  methods  (eg:  Statistical  Analysis  by 
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Cantwell,  i960). 

xhe  latter  method  was  adopted  mainly  because  the  data 
processing  could  be  done  entirely  by  electronic  means* 
oreover ,  this  method  does  not  require  the  lengthy  record 
that  is  required  for  (a),  nor  the  large  number  of  electronic 
circuits  that  is  required  for  (b). 

A  brier  description  of  the  order  of  data  processing  is 
as  follows.  First,  the  signal  (in  analog  form)  on  magnetic 
tape  was  converted  into  digital  form  by  means  of  the  play¬ 
back  circuit  shown  in  Figure  11.  This  digitized  data 
(punched  on  paper  tape)  was  then  fed  into  the  IBM  1620.  A 
program  was  written  to  calculate  first  the  power  spectra 
(proportional  to  volts^)  at  appropriate  frequencies,  and 
second  ,  the  apparent  resistivities  at  these  frequencies. 

Outline  of  Fathematical  Steps 

General 

Basically,  the  mathematical  operations  consist  of  two 
steps.  The  first  step  consists  in  determining  the  auto¬ 
covariance  C( T  )  at  lag t  ;  this  is  a  statistical  operation. 
The  next  step  consists  in  taking  the  Fourier  transform  of 
C(  nr  )  to  obtain  the  power  spectrum  P(f);  the  latter  is 
proportional  to  the  square  of  the  amplitude  at  frequency  f. 
Expressions  for  the  above  two  quantities  are  given  for  the 
case  of  a  continuous,  infinitely  long  record.  Modif ications 
to  these  expressions  are  then  made  when  the  record  is  of 


finite  duration  and  in  digital  form.  No  derivations  will 
be  given. 


Continuous,  Infinite  Record 

Expressions  for  C(T  )  and  P(f)  for  the  case  of  a 
continuous,  infinitely  long  record  will  be  given.  Suppose 
X(t)  is  the  voltage  between  two  points  on  the  earth’s 
surface.  The  autocovariance  at  lag  t  is  then  defined  as: 


C(  x)  =  Ove  [  X  ( t  +r)  .  X  (t  )  ] 


=  Jti 


m  zT  I  XCt+Xj.Xitjdt 
t  —  o o  _T 
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where  the  bar  denotes  complex  conjugacy.  In  this  case 
X(t)  is  real  and  the  bar  may  be  omitted.  When  X  s  0, 

C(0)  is  called  the  variance  and  is  given  by: 

T 

C(o)  =  / Lrn  —=  f  X(t )  .  X  ( t )  cl  t  . 

T-  oo  I  J  •  •  •  •  IV  -  d 
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Inspection  of  IV  -  2  shows  that,  if  the  resistance  between 
the  electrodes  be  taken  as  1  ohm,  then  C ( 0 )  is  a  measure  of 
power.  In  fact,  C ( 0 )  represents  the  average  amount  of 
power  dissipated  per  unit  time,  that  is  the  sum  of  the 
powers  of  all  the  component  frequencies.  Hence  C ( X  )  may  be 
written  as: 
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where  P(f)  is  referred  to  as  the  spectral  density  or 


simply  the  power  spectrum.  The  Fourier  transform  of 
IV  -  3  is: 


P<f)  =  f  C u) .  e 
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An  alternative  expression  for  p(f)  is,  using  IV  -  3 
and  IV  -  Ip : 
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The  power  spectrum  curve  with  suitable  modifications, 
usually  has  a  reasonably  simple  shape  which  represents 
fairly  accurately  the  actual  physical  situation  (Blackman 
and  Tukey,  page  88). 

When  the  unmodified  or  raw  power  spectrum  is  calculated 
from  IV  -  the  resulting  graph  has  a  very  irregular  fine 
structure  that  has  little  in  common  with  the  actual  physical 
situation  (H.  P.  F.  Swinnerton-Dyer,  (1962).  In  order  to 
avoid  this  undesirable  situation,  the  refined  power  spectrum 
can  be  calculated  by  using  IV  -  1  and  IV  -  ip  together  with  a 
"smoothing  function".  Hence  a  smoothed  power  spectrum  can 
be  calculated  directly. 

Finite  Record 

T.iOdlf ications  to  the  above  expressions  become  necessary 
when  the  records  are  taken  over  a  finite  time  interval,  Tn 
say.  C  (  TI  )  must  then  equal  zero  when  )  T  I  >  Ik  . 


A  function 


which  performs  this  operation  and  which  also  acts  as  a 
smoothing  function  is  the  following: 


IV  -  6 


o 


This  function  was  used  in  the  determination  of  P(f)  in 
this  thesis. 

Digitized  Record 

Digitization  of  the  data  introduces  an  effect  called 
"aliasing”.  The  effect  of  the  latter  is  to  limit  the 
upper  frequency  at  which  reliable  values  for  P(f)  are 
obtainable.  Suppose  readings  of  the  amplitude  of  a 
magneto  -telluric  signal  are  taken  at  uniformly  spaced 
intervals  t  *  0,  at,  ...nAt.  Then  aliasing  is  said  to 
occur  at  frequency  f^,  when  P(f-j_)  represents  not  just  the 
power  at  this  frequency,  but  in  addition,  the  sum  of  the 
powers  at  f-^  2  ,  f^  ±  JL  ,  etc.  In  order  to  ensure  that 

this  situation  does  not  occur,  an  upper  frequency  limit, 
fjj  say,  must  be  selected  so  that  all  the  desired  information 
is  contained  in  the  frequency  range  •  This  upper- 

limit  frequency  is  referred  to  as  the  "folding"  or  "Nyquist" 
frequency  and  is  given  by: 


IV  -  9 


where  is  the  time  interval  between  readings 
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Since  the  signal  response  of  the  amplifying  circuits 
was  poor  above  1  cps,  an  appropriate  value  for  a  t  would 
be  i  second,  so  that  f ^  =  1  cps.  Thus  P(f)  can,  in 
principle,  be  obtained  for  frequencies  less  than  1  cps. 

Once  a  value  for  is  selected,  then  the  power 
spectrum  can  be  calculated  at  equi-spaced  points  at 
frequencies : 


.  .IV  -  10 


where  m  is  the  number  of  lags.  The  larger  m  is,  the 

> 

greater  the  detail  that  appears  in  the  power  spectrum 
curve;  however,  this  detail  becomes  less  reliable.  An 
upper  limit  on  m  is  that  it  must  be  much  less  than  n,  the 
total  number  of  readings. 

Procedure  for  determining  power  spectrum  from  a 
digitized  record  of  uniform  spacing  of  finite  length  is 
as  follows: 

(1)  evaluate  the  mean  lagged  products  with  lag 


interval  AH  -  hat  by  using: 


IV  -  11 


(2)  evaluate  raw  power  spectrum  ( volts^/cps.2  toa  t  ) 
from: 


.IV  -  12 
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for  r  s  0,  1,  . .  .m 

(3)  evaluate  refined  power  spectrum  from: 

U  r  -  4"  Vr.,  +  I  Vr  +  Vr+  ,  .  .  ... IV  -  13 


Power  Spectrum  Graphs 

Graphs  of  power  spectrum  (volts^/cps  x  2m  4 t)  versus 
frequency  were  plotted  on  semi-log  paper  in  Figures  17  to 
21.  Each  figure  shows  the  power  for  the  E  component, 
represented  by  a  solid  line,  together  with  the  orthogonal 
H  component,  represented  by  a  dotted  line.  The  power 
spectra  curves  shown  are  for  the  amplified  signals. 

Properties  of  these  graphs  common  to  most  of  these 
curves  are  the  sharp  peak  in  the  0.01  to  0.03  cps 
frequency  range  and  the  steep  slopes  in  both  directions 
from  this  peak.  At  frequencies  above  0.2  cps,  most  of  these 
curves  show  a  very  jagged  fine  structure  which  is  the  result 
of  noise  (mainly  from  the  electronic  equipment  but  also 
from  the  signal).  The  power  in  this  frequency  range  is,  in 
general,  so  small  as  to  be  masked  completely  by  this  noise. 
The  noise  level  of  the  electric  component  (with  these 
amplifiers)  is  generally  greater  than  the  corresponding 
magnetic  noise  level  by  a  factor  of  10  (see  Figure  17b). 

Most  graphs  show  a  good  correlation  between  the 
electric  and  corresponding  orthogonal  magnetic  component. 

The  only  curves  which  show  an  increase  in  power  as 
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the  frequency  increases  from  0.2  cps  are  the  Bon  Accord 
curves.  A  second  maximum  point  in  those  curves  occur  at 
about  0.52  cps. 

The  curves  which  show  the  least  drop-off  in  power 
from  the  main  peak  at  0.02  cps  are  the  Cooking  Lake  power 
spectra.  These  curves  are  the  only  ones  which  yield 
reliable  information  in  the  frequency  range  0.2  to  0.5  cps. 

The  Irricana  power  spectra  shows  a  sharp  drop-off  in 
power  as  the  frequency  decreases  from  0.03  cps.  Consequently 
an  interpretation  at  great  depths  was  not  possible  because 
of  this  absence  of  power  at  long  periods. 

Since  the  maximum  frequency  at  which  useful  information 
is  obtainable  from  these  power  spectra  is  approximately  0.5 
cps  or  a  period  of  2  seconds,  the  shallowest  depth  at  which 
subsurface  structure  can  be  interpreted  from  these  curves 
is  2  kins  (skin  depth  p  &  ^  [iZpT  *»*)  •  This  value  can  be 
obtained  from  Table  1  when  the  surface  resistivity  is 
taken  as  10  ohm  meters  and  the  period,  2  seconds.  In 
central  Alberta,  the  resistivity  of  the  surface  layer  is 
just  under  10  ohm  meters. 
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CHAPTER  V 


RESULTS  AND  INTERPRETATION 

Calculation  of  Experimental  Resistivity  (in  practical  units) 

The  electric  field  intensity,  say  E__,  between  electrodes 

A 

with  a  separation  of  d(ft)  is  given  by: 


This  expression  gives  E  in  units  of  millivolts/kilometer. 

A 

The  amplification  is  obtainable  from  Figure  10(a);  the  power, 
P(f),  from  Figures  17  to  21. 

The  corresponding  magnetic  field,  II,  is  given  by: 


The  apparent  resistivity  pt  is  then  calculated  iron: 


(in  ohm  meters) 


Experimental  Apparent  Resistivity  Curves 

An  examination  of  Figures  23  to  2 (  show  the  lolloping 

features : 

(2.)  low  apparent  resistivities  in  the  short  period 
range  (2  to  10  seconds). 

(2)  sharp  rise  in  the  curves  when  the  period  increases 
from  10  to  100  seconds. 

(3)  a  drop  in  some  of  the  curves  in  the  region  100  to 


300  seconds. 

(4)  for  some  locations,  corresponding  pairs  of  curves 


(Ex>IIy  and  Ey,  Hx)  converge  as  the  period  becomes 
shorter. 

(5)  Tor  some  locations,  corresponding  pairs  of  curves 
tend  to  diverge  as  the  period  increases  from  10 
to  100  seconds. 

(6)  the  period  range  where  useful  information  is 
obtainable  from  a  curve  correlates  very  closely 
to  the  range  where  the  corresponding  power 
spectrum  curve  lies  above  the  noise  level. 

Five-layer  Case 

These  experimental  apparent  resistivity  curves  yield  no 
useful  information  for  periods  less  ^han  2  seconds;  this 
corresponds  to  a  skin  depth  of  2  kilometers  for  central 
Alberta.  Hence  an  interpretation  is  not  possible  for  depths 
shallower  than  this.  Moreover,  the  shape  of  the  apparent 
resistivity  curves  is  not  known  for  waves  whose  periods  are 
less  than  2  seconds  (that  is  whose  skin  depths  are  less  than 
2  kilometers).  In  Figure  22  are  shown  actual  earth  models 
together  with  corre sponding  apparent  resistivity  curves. 

The  data  for  these  models  were  obtained  from  electric  logs: 
the  locations  of  the  Leduc  and  Villingdon  wells  are  shown 
in  Figure  3«  The  resistivities  shown  for  each  layer  are 
average  values  estimated  by  a  visual  inspection  of  the  well 
logs;  no  attempt  was  made  to  obtain  the  true  bed  resistivity. 

These  particular  well  sites  were  chosen  because  their 
resistivities  for  the  upper  Cretaceous  are  extreme  cases; 


66 


the  maximum  resistivity  for  this  layer  (as  determined  from 
electrical  logs)  is  8  ohm  meters  and  the  minimum,  approximately 
[|_  ohm  meters.  The  resistivities  for  the  upper  Cretaceous  at 
the  5  recording  sites  lie  in  between  these  values.  In  addition, 
the  resistivities  for  the  upper  Devonian  at  Leduc  and 
VYillingdon  are  much  larger  than  those  at  the  recording  sites. 
These  high  resistivities  result  because  of  the  presence  of  the 
highly  resistive  reefs  in  the  Devonian  formation  at  these 
two  sites. 

Consequently,  one  might  expect  the  experimental  apparent 
resistivity  curves  to  lie  in  between  the  Leduc  and  Willingdon 
curves  in  the  period  range  0.01  to  1  seconds  and  to  lie  below 
the  Leduc  curve  in  the  period  range  1  to  10  seconds. 

From  this  magneto- telluric  study,  resistivities  for  the 
sedimentary  basin  were  obtained  which  agree  with  those  from 
electrical  logs.  The  fact  that  different  resistivities  for 
the  top  layer  (of  the  sedimentary  basin)  were  obtained  at 
different  locations  is  not  inconsistent  with  electrical  log 
values,  which  are  shown  in  Table  III. 

For  the  crystalline  basement,  the  resi stivi ties  were 
found  to  vary  from  600  uo  10^  ohm  meters.  However,  when 
the  value  10^  is  replaced  by  10^  in  the  theoretical  models, 
there  is  no  appreciable  difference  in  the  shape  of  the 
corresponding  apparent  resistivity  curve.  Hence  the  resist¬ 
ivity  for  the  precambrian  can  be  thought  of  as  varying  from 
600  ohm  meters  to  a  very  large  value,  say  10,000  ohm  meters. 
The  fact  that  the  resistivity  of  the  Precambrian  was  found 
to  vary  with  location  is  consistent  with  previous  geophysical 
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and  petrological  studies  oi'  the  basement.  Electrical  logs, 
from  wells  :n  central  Alberta  which  have  penetrated  into  the 
basement,  indicate  a  sharp  rise  in  resistivity  at  the  lower 
boundary  of  the  sedimentary  basin.  Moreover,  the  logs  indicate 
the  resistivity  to  vary  from  well  to  well  and  vary  from  100  to 
1,000  ohm  meters.  However,  these  are  superficial  values 
because  the  maximum  depth  of  well  penetration  into  the  base¬ 
ment  is  only  several  hundred  feet.  Hence  the  resistivity  of 
the  Precambrian  as  a  whole  can  be  expected  to  be  greater  than 
the  surface  values  because  the  moisture  content  of  the  rocks 
decreases  with  depth.  This  is  in  agreement  with  the  results 
of  this  investigation.  In  addition,  the  fact  that  correspond¬ 
ing  experimental  curves  (at  given  location)  tend  to  diverge 
is  suggestive  of  an  anisotropic  crystalline  basement. 

A  rather  significant  finding  is  that  no  electrical 
discontinuity  could  be  detected  at  the  depth  where  the 
Mohorovicic  discontinuity  is  expected  to  be.  This  suggests 
that  the  Mohorovicic  discontinuity,  which  is  a  seismic 
phenomena,  cannot  be  detected  by  the  magneto- telluric  method. 

An  electrical  discontinuity  was  detected  at  a  depth  of  00  km. 

Many  of  the  experimental  curves  tend  to  bend  over  at 
the  long  periods,  thus  indicating  the  presence  of  a  conducting 
substratum.  However,  other  factors  could  produce  the  same 
effect.  As  mentioned  in  Chapter  II,  the  effect  of  a  source 
of  finite  lateral  extent  is  to  decrease  the  apparent  resist¬ 
ivity  for  these  long  periods.  Hence  both  a  finite  source  as 
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well  as  a  highly  conducting  mantle  could  effect  a  bend-over  of 
the  experimental  curve  at  long  periods. 

An  interesting  feature  of  the  Cooking  Lake  experimental 
curves  is  that  one  curve  (Ex,  H  )  starts  to  bend  over  at  a 

J 

period  of  100  seconds,  whereas  the  other  (Ey,  Hx)  does  not. 

A  similar  situation  was  observed  on  the  experimental  curves 
obtained  at  Meanook  by  Niblett  and  Sayn-Wittgenstein  (i960). 
However,  their  (Ex,  H  )  curve  exhibits  this  feature  at  a 
much  longer  period  (close  to  1,000  seconds).  A  possible 
explanation  of  this  rather  unexpected  feature  could  be  a 
distribution  of  sources  in  such  a  way  as  to  modify  only  the 
(Ex,  Hy)  component  of  the  apparent  resistivity  at  long 
periods . 

Layer  of  Variable  Conductivity 

In  the  interpretations  so  far,  the  resistivity  for  the 
Precambrian  was  assumed  to  be  independent  of  depth.  However, 
the  Bon  Accord  (Ex,  Hy)  curve  (Figure  25)  gives  indication  of 
a  Precambrian  with  a  resistivity  which  increases  with  depth. 
Interpretation  by  Bossy  and  DeVuyst  (1959)  of  Belgium  of 
magneto- telluric  records  taken  at  Dourbes,  indicate  that  a 
composite  model  consisting  of  a  surface  layer  of  constant 
resistivity  overlying  a  layer  in  which  the  resistivity  varies 
continuously  with  depth  seems  to  represent  their  experimental 
results.  However,  since  only  one  curve  gives  indication  of 
this  feature,  no  interpretation  along  the  lines  of  Bossy  and 
DeVuyst  was  attempted. 
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TABLE  III 

APPARENT  RESISTIVITY  OP  TEE  SEDIMENTARY  FORMATIONS  OBTAINED 

PROM  ELECTRIC  LOGS 


Apparent  resistivity  in  ohrn  meters 


PERIOD 

FORMATION 

LEDUC 

mud  res.  2.8n.m. 
at  45°F 

ARDROSSAN 

1.2  mm 

at  60°F 

WILLINGD0N 

1.5^™ 

at  6i|°F 

Upper 

Edmonton 

8 

«» 

Cretaceous 

Lea  Park 

7 

6 

6 

Lower 

Viking 

3 

3 

3 

Cretaceous 

Joli  Fou 

7-10 

3 

3 

Nisku 

6o  -  300 

25 

- 

Upper 

Ireton 

10 

6 

7 

Devonian 

Reef 

Beaverhill  Lake 

50 

30  -  150 

15  -  ioo 

150 

10  -  30 

Middle 

Devonian 

10  -  30 

O 

u\ 

1 

o 

C\l 

25  -  50 

Upper  and 
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Middle 

10 

7-13 
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Pip.  23.  1  Experimental  and  theoretical  apparent  resistivity 
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LEGEND 

Solid  line  —  Experimental  Apparent  Resistivity  Curves 

A  and  o  —  Theoretical  Curves  and  cor  re  spond  i  ng 
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Pig.  25.  Experimental  and  theoretical  apparent  resistivity 
curves  together  with  models  for  the  theoretical 
curves . 
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CONCLUSION 

Theoretical  apparent  resistivity  curves  were  computed 
for  the  five-layer  case.  This  was  done  in  order  to  determine 
the  influence  the  sedimentary  basin,  which  consists  of  rocks 
from  five  geological  time  divisions,  has  on  the  shape  of 
the  apparent  resistivity  curves. 

A  convenient  program  has  been  devised  for  obtaining 
experimental  apparent  resistivity  curves.  The  data  is 
recorded  on  magnetic  tape.  Then  by  means  of  a  playback 
system,  this  data  in  analog  form,  is  converted  into  digital 
form.  This  digitized  data,  in  conjunction  with  programs 
designed  to  obtain  power  spectra  and  apparent  resistivity 
curves,  were  fed  into  the  electronic  computor  (IBM  1620  at 
the  University  of  Alberta) •  Thus  the  entire  data  processing 
was  done  by  electronic  means;  hence  the  most  time-consuming 
part  of  the  entire  operation  has  been  eliminated. 

Interpretation  of  the  experimental  apparent  resistivity 
curves  yield  values  for  the  resistivity  of  the  sedimentary 
basin  which  are  in  agreement  with  those  from  electrical  logs. 
Resistivities  ranging  from  6  to  9  ohm  meters  were  obtained 
for  the  near-surface  sedimentary  layers  and  values  from  6 
to  £0  ohm  meters  for  those  layers  near  the  Trecarnbrlan. 

The  results  of  this  investigation  indicate  the 
crystalline  basement  to  be  anisotropic.  However,  more 
recordings,  which  will  either  confirm  or  repudiate  this 
finding,  are  desired. 
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The  analysis  indicates  that  there  is  no  electrical 
discontinuity  at  the  Mohorovicic  discontinuity.  However, 
the  results  are  suggestive  of  an  electrical  discontinuity 
at  a  depth  of  80  km,  which  is  roughly  twice  the  depth  to 
the  Moho . 

At  the  long  periods  which  are  required  for  a  substratum 
investigation,  source  effects  may  be  significant.  Research 
into  this  aspect  is  presently  being  carried  out  at  the 
University  of  Alberta.  In  particular,  research  is  being 
conducted  to  determine  the  intrinsic  impedance  of  the  earth 
for  dipolar  fields  resulting  from  dipoles  (located  above  the 
earth’s  surface)  of  various  types  and  orientations. 
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APPENDIX 

Section  I 


PRACTICAL  UNITS 

H  measured  in  gammas (V) 

E  measure  d  in  m  , 

Length  measured  in  km 
p  measured  in  ohm  meters 


CONVERSION 

1  em  unit  =  10^  Y 
1  em  unit  =  |  "'/Km 
i  em  unit  r  10”^  km 
1  em  unit  r  10~^  ohm  meters 
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